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,__! ■ Abstract 

yo ■ The nature of the classical canonical phase-space variables for grav- 

O | ity suggests that the associated quantum field operators should obey 

2j ■ affine commutation relations rather than canonical commutation rela- 

— > . tions. Prior to the introduction of constraints, a primary kinematical 

rj-i representation is derived in the form of a reproducing kernel and its as- 

J ' sociated reproducing kernel Hilbert space. Constraints are introduced 

&0- following the projection operator method which involves no gauge fix- 

J> | ing, no complicated moduli space, nor any auxiliary fields. The result, 

which is only qualitatively sketched in the present paper, involves an- 
other reproducing kernel with which inner products are defined for the 
physical Hilbert space and which is obtained through a reduction of 
the original reproducing kernel. Several of the steps involved in this 
general analysis are illustrated by means of analogous steps applied to 
one-dimensional quantum mechanical models. These toy models help 
in motivating and understanding the analysis in the case of gravity. 

1 Introduction 

General relativity is, in certain ways, fundamentally different than most other 
physically relevant classical field theories, and the same remark applies to at- 
tempts to provide associated quantum formulations. The space-time metric 



S 



g^uix), x G IR 4 , //, v — 0,1,2,3, possesses a signature requirement that is 
incompatible with the space of metrics being a linear vector space. [] The 
inverse metric g a ^(x), defined so that g <T(J '(x)g flv (x) = 5°, is classically trivial 
but it is quantum mechanically challenged since the left-hand side involves 
the product of two operator-valued distributions. Moreover, the spatial and 
temporal constraints that hold at each space-time point classically close al- 
gebraically, but they exhibit an anomaly when quantized. In effect, this fact 
changes the constraints from first class (classically) to second class (quantum 
mechanically). And, of course, there is the well-known fact that unlike other 
theories which take place on a fixed space-time stage, the theory of gravity 
involves the dynamics of the space-time stage itself. The purpose of this 
article is to discuss some basic issues surrounding quantum gravity from a 
viewpoint different than traditional ones.Q 

Let us outline the general approach we shall adopt. First, we focus on 
basic kinematics and the quantum theory of positive definite, 3x3 matrix- 
valued field variables and associated noncanonical "conjugate" field variables, 
designed to offer an initial class of coherent states and coherent-state in- 
duced Hilbert space representations. In this step it is noteworthy that, be- 
sides the signature issue for the metric, the existence of an operator g^ k (x), 
j, k = 1,2, 3, inverse to gu{x) is shown, such that, when suitably defined, the 
equation g jk (x)gki(x) = 5{ is fulfilled. Second, we introduce the spatial and 
temporal constraints in the projection operator approach recently developed 
by the author and others H, ||, ||, 0. This procedure has the advantage of 
working entirely with the classical degrees of freedom, including the c-number 
Lagrange multipliers — specifically, the lapse and shift functions. It is not nec- 
essary to introduce additional fields (e.g., ghosts with false statistics), nor 
choose gauges, nor pass to moduli spaces, etc. Initially, the constraints are 
imposed in a regularized fashion. Subsequently, the removal of the regu- 
larization is analyzed, a process which often involves an automatic change 
of Hilbert-space representation. Assuming that the limit removing the reg- 
ularization exists, the physical Hilbert space that arises is then, generally 
speaking, best described as a reproducing kernel Hilbert space that emerges 

1 Although, we assume a 3 + 1 theory of gravity for illustrative purposes, it is straight- 
forward to generalize to an s + 1 theory as well, s > 1. 

2 The closest work in spirit to that discussed in this paper is that of the author 0], 
Isham and Kakas Q, and especially Pilati |!|. See Sec. 5 for an extensive discussion. We 
do not directly comment on current schemes for quantizing gravity. 



from the reduction of the original reproducing kernel. 

Before we undertake any discussion of gravity, however, we sketch in Sec. 2 
the key concepts as applied to some simple, few degree-of-freedom systems. 
In Sec. 3 we construct a suitable kinematical framework for quantum gravity, 
while in Sec. 4 we analyze the introduction of constraints. Section 5 contains a 
general discussion about operator representations and constraints in relation 
to reparameterization invariance. In Part II of this work, the analysis of 
gravitational constraints is discussed in detail. In addition, the classical 
limit of the affine gravitational quantum theory developed in Sees. 3-4 will 
be discussed and compared with classical gravity. 

2 Elementary Illustration of Key Concepts 

As is generally well known, there is only one irreducible representation up to 
unitary equivalence of canonical, self-adjoint operators P and Q satisfying 
the Weyl form of the canonical commutation relations. This representation, 
equivalent to the Schrodinger representation, implies that the spectrum of 
both P and Q cover the whole real line. Such operator degrees of freedom are 
appropriate for many systems with a finite or an infinite number of degrees 
of freedom, but they are inappropriate for gravity. The reason for this is 
that the classical 3x3 metric is strictly positive definite and the associated 
quantum field operator cannot be represented by an operator whose spectrum 
is unbounded above and below. Instead, of the usual relation [Q, P] = i, with 
h = 1, one is led to consider an affine commutation relation ||, which for a 
single degree of freedom takes the form 

[Q,D] = iQ. (1) 

Here D = (PQ + QP)/2 denotes the dilation operator, and it follows || |I0| 
that solutions of the affine commutation relations exist with irreducible, self- 
adjoint operators D and Q for which — and this is the important part — Q > 0. 
(There are two other inequivalent self-adjoint solutions, one where Q < 0, 
which is rather like the representation of interest, and another for which 
Q = II . Neither of these representations will be of interest in this article.) 
Even though the operator P is only a symmetric operator which has no self- 
adjoint extension, the introduction of the self-adjoint operator D provides 
the substitute commutation relation given above. These two commutation 



relations are not in conflict since the affine commutation relation follows 
directly from the Heisenberg commutation relation simply by multiplication 
of the latter by Q. We note that an analog of the affine variables will be used 
in the case of the gravitational field to maintain the positivity of the local 
quantum field operator for the 3x3 spatial metric. 

Continuing with the one-dimensional example, and based on self-adjoint 
operators that satisfy the affine commutation relation, let us introduce affine 
coherent states [|ll|], \p,q) G $), defined by the expression 



\p, q) = e ipQ e~ iXn{q)D \rf) , -oo < p < oo , < q < oo . (2) 

Here, the fiducial vector |?7) is chosen to satisfy several conditions, which, 
using the shorthand ((•)) = (77K • )\r]), are specifically given by 

(Q- 1 )=C<oo, (1) = 1, (Q) = l, (D) = 0, (P) = 0. (3) 

The first condition is required, while the remaining conditions are chosen 
for convenience. The coherent states also admit a resolution of unity |iT 
expressed in the form 

1 = / \P, <l) (P, q\ dr(p, q) , dr{p, q) = dp dq/2nC , (4) 

integrated over the half plane R x 1R + . 

In particular, diagonalizing the self-adjoint operator Q = f™x\x)(x\ dx 
in terms of standard Dirac-normalized eigenvectors, leads to a representation 
for the coherent-state overlap given by 

(p,q\r,s) = ( v \e iln ^ D e- ipQ e irQ e- iln ^ D \ V ) 
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= (qs) ' rj(x/qye-' x(p - r) 7 1 (x/s)dx } (5) 

Jo 

where the fiducial function rj(x) = (x\r)) denotes the Schrodinger representa- 
tion of the fiducial vector |?7). It is important to observe, for some suitable 
function F, that 

(p,q\r, s) = F(q,p-r, s) , (6) 

namely, that p and r universally enter in the form p — r. It is also clear that 
(p,q\r,s) defines a continuous, positive- definite function, which, apart from 



the continuity, means that 

N 

22 a *n a m (Pn, q n \Pm, <?m) > (7) 

n,m=l 

for arbitrary complex {a n }^ =1 and real {p n , q n }n=i sequences, with N < oo. 
The function (5) may be taken as the reproducing kernel for a reproducing 
kernel Hilbert space [I2|. Note that the information in (p,q\r,s) is enough 



to recover r/(x) apart from an overall constant phase factor. Thus different 
fiducial functions (not related by a constant phase factor) generate distinct 
reproducing kernels. Since each reproducing kernel Hilbert space has one 
and only one reproducing kernel [0], it follows for different i](x) that the 



Hilbert space functional realizations are completely disjoint except for the 
zero element. Basic elements of a dense set of vectors in each such Hilbert 
space are given by continuous functions of the form 
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V>(P> Q) = ^2 a n(p, q\p n , q n ) (8) 
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defined for arbitrary complex {a n }^ =1 , and real {p n , q n }n=i sequences, with 
N < oo. Let a second such function be given by 

j 

<f>(p, q) = Yl & & q % ' 3* ) ( 9 ) 

defined for arbitrary complex {[3j}j =l and real {Pj,qj}j = i sequences, with 
J < oo. The inner product of two such vectors is then defined |12[] to be 
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(</#> = (V(-,-) ,0( v )) =Y,Y,<l 5 ^qn\p j ,q j ) , (10) 

n=l j=l 

which when \<f>) = \ip) is, by definition, nonnegative. The resultant pre- 
Hilbert space is completed to a (reproducing kernel) Hilbert space C by in- 
cluding all Cauchy sequences in the norm |||'0}|| = Vv^lW as N tends to 
infinity. Lastly, we note that the space of functions appropriate to one repro- 
ducing kernel is identical to the space of functions appropriate to a second 



reproducing kernel that is just a constant multiple of the first reproducing 
kernel. This fact does not contradict the uniqueness of the reproducing kernel 
for each Hilbert space because strictly different inner products are assigned 
in the two cases. Of course, the foregoing discussion applies quite generally 
and is not limited to any one sort of reproducing kernel. 

When the states \p, q) form a set of coherent states — as we assume in the 
present case — the inner product has an alternative representation given by a 
local integral of the form 

(V#) = / i>{p,q)*<f>{p,q)dT(p,q) , (n) 

expressed in terms of ^{p.q) = (p,q\4>) and (fi(p,q) = (p,q\4>)- It follows 
that this formula holds for all elements of the completed Hilbert space C, 
and, moreover, every element of the so-completed space is a bounded and 
continuous function, the collection of which forms a rather special closed 
subspace of L 2 (R 2 , dr). 

Thanks to the coherent-state resolution of unity, it follows that the co- 
herent state overlap function satisfies the integral equation 

(p",q"\p',q') = J(p",q"\p,q)(p,q\p',q')dr(p,q) , (12) 

a basic relation, which, if it was established as a first step for the continuous 
function (p",q"\p',q') [= (p',q'\p",q")*}, guarantees the existence of a local 
integral representation for the inner product of two arbitrary elements in the 
associated reproducing kernel Hilbert space. Several useful properties follow 
from this reproducing property. For example, repeated use of the resolution 
of unity leads to the fact that 

/« Li Li 

•• / Y[(Pi+i, Qi+i\Ph Qi) IJ rfr (^'^) ' ( 13 ) 

J 1=0 1=1 

in which we have identified p", q" = pL+i,qi+i and p',q' = po,qo- In turn, 
making an (unjustified!) interchange of the (continuum) limit with the inte- 
grations, and writing for the integrand the form it would assume for contin- 
uous and differentiable paths, gives rise to the suggestive but strictly formal 
expression JT3 1 

(p'VbW) = je^I T o^^ dt Vr(p,q) , (14) 



which determines a formal path integral representation for the kinematics 
that applies for any T > 0. Thus, the existence of a coherent state resolution 
of unity is the necessary condition to introduce a traditional coherent state 
phase-space path integral representation for the kinematics, which specifically 
leads to the reproducing kernel. A path integral for the kinematics is a 
necessary prerequisite to obtain a path integral for the dynamics. While less 
than ideal, the formal path integral itself may be used as a starting point 
for quantization. Even though the formal nature of the path integral renders 
it basically undefined, one may always (re)introduce a regularization by a 
lattice-limit formulation (as above), using suitable ingenuity to choose an 
acceptable integrand. This procedure is more or less standard by now. 

However, it must be appreciated that reproducing kernels for reproducing 
kernel Hilbert spaces are not required to fulfill a (positive) local integral 
representation for the inner product. In cases where the integral for the 
resolution of unity does not exist, one must accept the inner product that 
is given directly by the reproducing kernel, which means that the integral 
relation (12) is replaced by 

((■r\p",q"} 3 (-\p / ,q')) = (p",q"\p / ,q / ). (15) 

When this is the case we say that {\p, q)} forms a set of weak coherent states 



14]], i.e., the elements of {\p, q)} span the Hilbert space fj, but do not ad- 
mit a local integral representation for the inner product of elements in the 
associated reproducing kernel Hilbert space. 

A simple example of a reproducing kernel Hilbert space without a local 
integral representation for the inner product is determined, for u", u' G R, 
by the reproducing kernel (u"\u') = exp[— (u" — u') 2 ]; here, one must use 
(( ' \ u ") 3 ( ' \ u ')) — ( u "\ u ')-> which is then extended by linearity and continuity 
to all Hilbert space vectors. 

A more relevant set of examples is given by the following discussion ap- 
plied to our simple model. Let a > —1/2, and choose rj(x) = N x a exp(— f3x). 
Here the factor N is fixed by requiring J \rj(x)\ 2 dx = 1. The two condi- 
tions (Q) = 1 and (Q^ 1 ) = C < oo, lead to (3 — | = a > 0; in this case 
C — 1 — 1/(2/?). In turn, the reproducing kernel is given explicitly J15) by 

r {qs)-^ 120 

\P,q\r,s) [ij g - 1 + s - 1)+i i (3 - 1{p _ r 



= exp(-2p\n{H( q -i + s-i) + i\^(p - r)]/(qs)-^}) ; 

(16) 

the second form is given for comparison purposes to the gravitational case. 
As long as (3 > \ it follows that the states \p, q) are a set of coherent states 
with a proper resolution of unity and therefore a local integral representation 
for the inner product exists. In this case, path integrals exist as lattice 
limits, and the whole situation seems familiar. On the other hand, if < 
(3 < |, the overlap function (p, q\r,s) defined above is still a positive-definite 
function and, therefore, it is a valid reproducing kernel which leads to an 
associated reproducing kernel Hilbert space; however, such a Hilbert space 
does not admit a local integral representation for the inner product in terms 
of the given representatives. Therefore, there is no conventional coherent 
state path integral for the kinematics, and thus also for the dynamics, in a 
reproducing kernel Hilbert space representation when < j3 < |. This lack of 
a conventional coherent state path integral representation may appear to be 
detrimental to any program to introduce quantization, dynamics, etc. — but 
there is hope. 

There is another way to generate the reproducing kernel for the given 
family of fiducial vectors which may be applied for all /3 > 0. Let us first 
focus on (5 > |. In that case, observe, by construction and using d p = d/dp, 
etc., that for every \tp) (Eft, 

B^( P ,q) = {-iq-% + 1 + /T 1 ? W(p,<?) = ° » ( 17 ) 

an equation which represents a (complex) polarization [|l(| of L 2 (R 2 ,g?t). It 
follows that the second-order differential operator A = \(3B^ B > 0, and 
therefore A can be used to generate a semigroup. In particular, for any 
T > and as v — >• oo, the expression e~ uTA becomes a projection operator 
onto the subspace C of solutions to the polarization equation [[[]]] . In a two 
degree of freedom Schrodinger representation — symbolized by \p,q), where 
(p, q) e R x R + and (p, q\r, s) = 5{p — r)5(q — s) — it follows, from a two- 
variable Feynman-Kac-Stratonovich path integral formula [ITH], that 



(p",q"\p',q')= lim (p",q"\e uTA \p',q') 



iuii M ! e^I o q ^P^ dt ~^^ 2 ^I o^' 1 ^ 2 ^^^^' 2 ^ 2 ^ 1 VpVq 
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= lim e» T ' 2 [ e-'fo*® dp{t) dW v {p, q) , (18) 

where W u denotes a two-dimensional Wiener measure with diffusion constant 
v, pinned at t — to p', q' and at time t — T to p", q", which is supported 
on a space of constant negative curvature R = —2/(3. It is noteworthy in 
(18) that the variable p(t) enters only in the form p(t), a fact which leads 
to the result depending only on the difference, p" — p'. For every v < oo, 
and with probability one, all Wiener paths in the given path integral are 
continuous, and, for purposes of coordinate transformations, it is convenient 
to adopt the (midpoint) Stratonovich rule to define the stochastic integral 
—fq(t) dp(t) since, in that case, the rules of the ordinary calculus hold. Such 
a representation is said to involve a continuous-time regularization flT9[| . 

Now we consider the case where < (3 < |. The solutions to (17) are, up 
to a factor, analytic functions, but they are no longer square integrable, as 
is clear from the fact that the large q behavior of (16) is controlled only by 
the factor q~@ . As a consequence, the operator A > has only a continuous 
spectrum. The family of operators e~ uTA is still a semigroup, and the ex- 
pression (p",q'\e~ uTA \p',q') still has the formal path integral representation 
given in the middle line of (18). However, as v — > oo, e~ uTA does not lead to 
a projection operator; instead we need to extract the germ of that semigroup 
as v —>■ oo. If we let E > denote continuum eigenvalues for the operator 
^B^B, we can write 

{p'\q"\e-» TA \p\q') = I \p" ,q"\E,v) e~^ E (E , v\p' , q') p(E,v) dE dv , (19) 



for some density of states p(E, v). Here, the variable v labels degeneracy for 
^B^B. For the sake of illustration, let us assume that p(E,v) ~ CE w ~p(v), 
w > — 1, for E <C 1. We base this assumption on the fact that there is no 
reason for an E'-dependent degeneracy for very tiny E. Thus, we consider 
the expression J{y) = [y (3T) W+1 / 'C T [ (w + 1), and are led to the fact that 

lim J{y) (p", q"\e~ uT \p', q') 

v— »oo 

lim j{p\ q"\E, v) e~ uTf3E (E, v\p', q') J{v) p{E, v) dE 

iy->ooj 

(p",q"\0,v)(0,v\p',q')p(v)dv. (20) 



In effect, this procedure has enabled us to pass to the germ of the semigroup. 
Observe that the rescaling factor is independent of the coherent state labels, 
and thus we are only making a //-dependent rescaling before the limit v — > oo 
is taken. If necessary, we can rescale our expression by letting J[y) — > 
J(u) = M(p", q") M(p', q') J{v) to achieve normalization without effecting 
its positive-definite character. In summary, for < (3 < |, we claim that 
instead of (18) we can write 



(p\q"\p\q')= lim J{ v ) (p", q"\e^ TA \p' : q' 



= lim 77 ! e ^^ o q ^P^ dt -^/ 2 ^I o^ lq ^ 2 P^ 2+ ^ q ^ 2( i^ dt VpVq 
u ->°° J 

== lim 7(i/) e uT/2 [ e _i /^ (t) dp(t) dW u (p, q) , (21) 



v— >oo 



Convergence in this case is initially regarded in the sense of distributions. To 
lead to the desired result, we appeal to analyticity (up to a specific factor) 
of the result in q~ l + i(3~ 1 p, analyticity (up to another factor) in the vari- 
able s" 1 — if3~ 1 r ) and dependence on p — r. Note that J(v) can always be 
determined self consistently by insisting that (p,q\p,q) = 1 for all (p,q). In 
simpler terms, we can always regard J(u) as part of the needed normalization 
coded into AT in the formal path integral expression. 

As will become evident later, various features of reproducing kernels, 
reproducing kernel Hilbert spaces, and associated rules for defining inner 
products illustrated above will carry over into the quantum gravity case as 
well. 

2.1 Operators and symbols 

In addition to the properties of the reproducing kernel Hilbert space, certain 
symbols associated with operators are important. Let us introduce the up- 
per symbol H(p, q) associated to the operator H,(P, Q) and defined, modulo 
suitable domain conditions, by the expression 

H(p,q) = (p,q\H(P,Q)\p,q) = (H(p + P/q,qQ)) . (22) 

For example, if 7i(P, Q) = P 2 — Q^ 1 denotes a quantum Hamiltonian, then 
H(p, q) = p 2 + (P 2 ) /q 2 -C /q. Since C = 0(1) (e.g., for large (3) and (Q) = 1, 
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it follows that (P 2 ) = 0(h 2 ). Observe that H basically agrees with the ex- 
pected classical Hamiltonian in the limit that ft — > 0, but prior to that limit 
H includes a quantum induced barrier to singularities in solutions of the 
usual classical equations of motion. We adopt the expression H(p, q) as the 
(h- augmented) classical Hamiltonian and refer to the connection between the 
quantum generator TC and the classical generator H as the weak correspon- 
dence principle |20| . In this way, the classical and quantum theories may 
both coexist, as they do in Nature. 

There is also another set of symbols that are important. We introduce the 
lower symbol h(p, q) which is related to the operator TC(P, Q) by the relation 

H(P, Q) = J h(p, q) \p, q) (p, q\ dr(p, q) . (23) 

For the one-parameter class of fiducial vectors leading to (16), it follows that 
a dense set of operators admit such a symbol for a reasonable set of functions. 
In the quantum gravity case, there are twin goals: (i) to ensure that 
the field operators of interest are well defined and locally self adjoint in the 
given field operator representation; and (ii) to choose locally self adjoint con- 
straint operators that have a weak correspondence principle which connects 
them with the desired form of the classical constraint generators (possibly h 
augmented) . 

2.2 Imposition of constraints 

We adopt the projection operator approach to the quantization of systems 
with constraints ||, [5|, |, 0. Let {$ a (P, Q)}^ =1 , A < oo, denote a set of 
constraints each given by a self-adjoint operator. Further assume that $•$ = 
£^ =1 ($ a ) 2 is also self adjoint. We define the (provisional) physical Hilbert 
space .fiphys = Ef), where E = E^ = E 2 is a uniquely defined projection 
operator, and in turn choose 

E = E($ • $ < 5(h) 2 ) , (24) 

where 5(h) is not a (^-function but a small, positive, possibly h dependent, 
regularization parameter for the set of constraints. As shown below, 5(h) is 
chosen so that E is the desired projection operator. This choice may entail 
a specific representation of the Hilbert space and a suitable limit as 5 — > 
to extract the germ of the projection operator. 

11 



We may illustrate this latter situation for the constraint $ = Q — 1, 
assuming initially that 5 < 1. In this case 

rl+S r- 

ME((Q - l) 2 < 5 2 ) |0) = y rfx y dtr(y) ^(x, y)*<t>{x, y) , (25) 

where a accounts for any degeneracy that may be present. When restricted 
to functions ip and <p in the dense set D, where 

D = { polynomial (x, y) e~ x2 ~ y2 } , (26) 

and rescaled by a suitable factor, the projection operator matrix elements 
lead to the expression 

(25)- 1 (^ o |E((Q-l) 2 <5 2 )|0 o ) 

/l+<5 /■ 

dx J da{y)^ {x,y)*(t> {x,y). (27) 

Now, as 5 — ► 0, this expression becomes 

^ (1, y)> (l, y) da(y) = ((^ , O )) . (28) 

Interpreting this final expression as a sequilinear form, one completes the 
desired Hilbert space by adding all Cauchy sequences in the associated norm 
|| \ip )) || = y/((i/j , i'of) ■ The result is the true physical Hilbert space in which 
the constraint Q — 1 = is fulfilled, and this example illustrates how con- 
straints are to be treated when $ ■ $ has its zero in the continuous spectrum. 
A second example of an imposition of constraints is given by $i = Q — 1, 
as before, along with $ 2 = D. This situation corresponds to second class 
constraints, and serves as a simple qualitative model of what occurs in the 
gravitational case. In this case 

E = E(D 2 + (Q-1) 2 <5(£) 2 ). (29) 

Here, the left-hand side of the argument can be regarded as a "Hamiltonian" , 
and the ground state |0') for such a system (nondegenerate, in the present 
example) can be sought. In particular, there are two positive parameters, 5' 
and 5", such that, for all S with 5' < 5 < 5", then 

E = E(Z) 2 + (Q- l) 2 < 5(h) 2 ) = |0'> (0'| . (30) 

12 



This is the desired choice to make for IE in the case where $ ■ $ has a discrete 
spectrum near zero that does not include zero. 

For completeness, if $ • $ has a discrete spectrum including zero, then it 
suffices to choose 

E = E($-$ = 0) = E($-$<5(ft) 2 ) , (31) 

where in the present case 5(h) > is chosen small enough to include only 
the subspace for which $ • $ = 0. This simple case does not seem to arise in 
quantum gravity. See (|, |7j for a discussion of gauge invariance. 

When the number of constraints is infinite, A = oo, as will be the case 
for a field theory, then a slightly different approach is appropriate. One form 
this takes is dealt with in Sec. 4. 



2.3 Appearance of time 

In a reparametrization invariant problem in quantum mechanics it is typical 
that dynamics is cast in the guise of kinematics at the expense of introducing 



an additional degree of freedom plus a first-class constraint; see, e.g., |H 
Let the resultant kinematical reproducing kernel with the extra degree of 
freedom be given by (p",q",s",t"\E\p',q', s',t'), where E is the projection 
operator enforcing the constraint. Next, reduce this expression, for example, 
as in the procedure 

(p", q", t"\p', q', t>) ee ff(p", q", s" , t"|E|p', q', s ' , t>) ds" ds> . (32) 

The result is a new positive-definite function that can be used to define a 
reproducing kernel Hilbert space. However, it may well happen that the 
following identity holds 

(p",q",t"\p>,q>,t>) = ({;;Y,q",t") 3 {;;-\p',q',t'}) 

= ((;;t\p",q",t") 3 (;;t\p',q',t')). (33) 

This equation means that the space spanned by the states \p, q, t) by vary- 
ing p, q and t is the same space spanned, by the same states, by varying 
p and q but with t held fixed at some value (e.g., t = 0). This situation 
implies that the states \p, q, t) are extended coherent states in the sense of 
imj , and, in particular, thanks to using canonical group coordinates in the 

13 



coherent-state parameterization, that \p, q, t) = exp(—iHt)\p,q,0) for some 
self-adjoint "Hamiltonian" 7i. The parameter t is then recognized as the 
"time" t. For an explicit example of how this procedure works in detail see 
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It is expected that a suitable time parameter will emerge in the extension 
of these ideas to the gravitational case. 

2.4 Matrix generalization 

Our preceding analysis has been confined to a single p and q and the asso- 
ciated affine quantum operators. In any generalization to the gravitational 
case, it will first be necessary to generalize the preceding discussion to 3 x 3 
(or more generally to s x s) matrix degrees of freedom and repeat an analysis 
similar to that of the present section. We do not include this discussion here 
since that is the subject of a separate work ||24|| . It is safe to say, apart from 
some technical details, that there are no special surprises in this generaliza- 
tion, and the basic concepts that we shall need are already present in the 
simplest case on which we have concentrated. 

3 Gravitational Kinematics 

3.1 Preliminaries 

Let us start with the introduction of a three-dimensional topological space S 
which locally is isomorphic to a subset of R 3 . Locally, we generally use three 
"spatial" coordinates, say x J , j = 1, 2, 3, to label a point in S. This labeling is 
nonsingular and thus one-to-one. Whether a single coordinate chart covers 
S depends on the global topological structure of S. Let us fix this global 
topological structure from the outset — for example, topologically equivalent 
to R 3 , S 3 , T 3 , etc. The theory of quantum gravity developed here does not 
engender topological changes of the underlying topological space S. Note 
this lack of topological change applies only to the space S. It is unrelated 
to any presumed "space" and/or "topology" associated with any quantum 
metric tensor, which, after all, is typically distributional in character. 

If the space S is such that more than one coordinate patch is required we 
arrange for the necessary matching conditions and rename the coordinates 
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within each patch by x\ j = 1, 2, 3, for some domain. We can also consider 
alternative coordinates, say ~x\ j = 1,2,3, which are also nonsingular. We 
admit only differentiable coordinate transformations such that the Jacobian 
[dx/dx] = &et(dxi /dx k ) ^ everywhere. The group composed of such 
invertible coordinate transformations is the diffeomorphism group. 

We can also introduce functions on the space S which in coordinate form 
may be denoted by f(x). A scalar function is one for which f(x) = f(x), 
while a scalar density of weight one satisfies b(x) = [dx/dx] b(x), or stated 
as a volume form, dV = b(x) d?x = b(x) d 3 x. Observe that it is not neces- 
sary to have a metric in order to have a volume form. Whether S is com- 
pact or noncompact, we assume that < b{x) < oo for all x and therefore 
< 6(x) _1 < oo for all x as well. These properties are still valid after a 
nonsingular coordinate change. Integrals of a scalar function take the form 
J f dV = j f(x) b(x) d 3 x and are invariant under any coordinate transforma- 
tion in the diffeomorphism group. 

In the ADM (Arnowitt, Deser, Misner) p5| canonical formulation of clas- 
sical gravity there are two fundamental fields gki(x) [= gik(x)] and n kl (x) [= 
n lk (x)]. The metric gki{x) transforms as a (two-valent covariant) tensor, while 
the momentum 7i kl (x) transforms as a (two-valent contravariant) tensor den- 
sity of weight one. Thus f gki(x) n kl (x) d 3 x [or even J gki(x, t) n M (x, t) d 3 x dt] 
is an invariant under diffeomorphism group transformations (on the spatial 
hyperspace, of course). Note well: The latter example pertains to a gener- 
alization of the former one including an additional independent variable t 
which possibly could be identified with (coordinate) "time". 

A metric is not an arbitrary tensor but is restricted to be positive definite. 
Specifically, for any real a? , j = 1,2,3, where S^ =1 (a J ) 2 > 0, it follows 
that a k gki(x) a 1 > for all x. As a consequence, the positive-definite (two- 
valent contravariant) tensor g kl (x) exists at each point and is defined so that 
g kl (x) gi m (x) = <5m- In addition, \/g(x) = y/det[gki(x)] > transforms as a 
scalar density of weight one. Thus, as is well known, \J 'g(x) d 3 x characterizes 
a volume form, but this choice ties the volume form to a specific metric, or at 
least to a specific class of metrics. This close association to specific metrics 
is something we would like to avoid, and it leads us to choose b(x) d 3 x as 
the preferred volume form. Of course, if b(x) = \Jg(x) everywhere in any 
coordinate system, then the volume form b(x) d 3 x is identical to the one based 
on a metric space and given by \/g(x) d 3 x. 
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3.2 Reproducing kernel — original Hilbert space 

A study of canonical quantum gravity begins with the introduction of metric 
and momentum local quantum field operators, which we denote by a k i(x) [= 
<Jik(x)] and fi kl (x) [= fi lk (x)}, respectively. For such fields one postulates the 
canonical commutation relations 

[(Tki(x), fi ra (y)] =iS r k a l 8(x,y) , 

[a ki (x), cr rs (y)] = , 

[//<(x),//%)] = 0, (34) 

with Sli = (8 k Sf + 5[<5|)/2. Since the right-hand side of the first equation 
is a tensor density of weight one, it is consistent that we define a k i(x) to 
be a tensor and fi rs (x) to be a tensor density of weight one. However, just 
as its one-dimensional counterpart, there are no local self-adjoint field and 
momentum operators that satisfy the canonical commutation relations as 
well as the requirement that {o~ k i(x)} > 0. To arrive at a suitable substitute 
set of commutation relations, we introduce, along with the local metric field 
operator a k i(x), the local "scale" field operator K k (x) which together obey 
the affine commutation relations [[[], ^ Q 

KOg). K t(y)] = i Wi 40) - S l K k(x)} <*fo y) , 

[o-m(x), K r s (y)} = i i [81 a ls (x) + 8\ a ks (x)} 5(x, y) , 

[a k i(x), cr rs (y)] = . (35) 

In these relations, o~ k i(x) remains a tensor, while K r s {x) is a tensor density of 
weight one under coordinate transformations. The local operators K r s {x) are 
generators of the GL(3,R)°° group [0, §], while the local operators a k i(x) are 
commuting "translations" coupled with the GL(3,R)°° group by a semi-direct 
product. The given affine commutation relations are the natural generaliza- 
tion of the one-dimensional affine commutation relation presented in (1). In 
the case of one degree of freedom, the affine commutation relations follow 
from the canonical ones, while in the case of fields this is, strictly speaking, 
incorrect. It is true that 

K l( x ) = Wki{x) fi lr (x) + fi rl (x)ai k (x)} R , (36) 

where the subscript R denotes an infinite multiplicative renormalized prod- 
uct to be defined later. However, the presence of an infinite rescaling means 
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that either the canonical or the affine set of commutation relations can hold, 
but not both at the same time. Since it is the affine commutation relations 
that are consistent with local self-adjoint operator solutions enjoying metric 
positivity, we shall adopt the noncanonical affine commutation relations. The 
choice of the affine commutation relations means that the canonical commu- 
tation relations do not hold, and therefore we are dealing with a noncanonical 
quantization of the gravitational field. 

Accepting the affine field operators as generators, we introduce a primary 
set of normalized affine coherent states each of which — in a deliberate abuse 
of notation — is defined by 

1 71- q\ = e *J ^'(^^kd^d,^ e -iJ^(x)K a r (x)d 3 x I \ /^7) 

for a suitable fiducial vector \rj) characterized below. In \n,g) the argument 
"tt" denotes the momentum matrix field n ab while "o" denotes the metric 
matrix field g a b- By all rights, the states in question should have been called 
\ir, 7}, but as we shall see, the overlap of two such states, for the featured 
choice of \rj) [cf. (70)], depends only on the matrix (for each point x) g = 
exp(7 T /2) exp(7/2) = {g a b}, where T means "transpose" .0 If the space S is 
noncompact, then, as smooth c- number fields, both n and 7 should go to zero 
sufficiently fast so that the indicated smeared field operators are indeed self- 
adjoint operators and generate unitary transformations as required. On the 
other hand, as we shall shortly see, this asymptotic behavior can, effectively, 
be significantly relaxed. 

The overlap of two such coherent states leads to an expression of the form 

(n",g"\n',g') = F(g",n"-n',g') (38) 

for some continuous functional F which depends only on the difference of 
the fields, n"(x) — vr'(x), an analog of which already occurred for the one- 
dimensional example. Whatever choice is made for the fiducial vector |n), 
the coherent state overlap function defines a continuous, positive-definite 
functional which, therefore, defines a reproducing kernel and its associated 
(separable) reproducing kernel Hilbert space C. By construction, therefore, 

3 Observe by this parameterization that g > as opposed to a traditional triad for which 
g > 0. To see that this distinction may possibly make a real difference see p5[ . We remark 
that the nonsymmetric matrix exp(7/2) would have relevance for spinor fields. 
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the set of coherent states {\tt, g)} span the Hilbert space $). As such they 
form a basis (overcomplete to be sure!) for S). Based on arguments to follow, 
we are led to the proposal [cf., (16)] that 

(n",g'V,g'} = 

(-0 fh( W3 1 f det iW kl ( x ) +g' k K*)} +i\b(x)- l [n" k \x) - 7T fkl (x)]} ^ 
6XP l z J^ x ) dx[n { (det[^(x)])V2(det[^(a;)]) 1 /2 )) 

(39) 

This equation is central to our analysis of quantum gravity. 

The coherent-state overlap (39) may be read in two qualitatively dif- 
ferent ways. Although arrived at on the basis that 7i"(x), tt'(x) — > and 
7" '(x) , y (x) — * as \x\ — »■ 00, the given expression exists for a far wider 
limiting behavior. In particular, suppose there is a fixed asymptotic behavior 
such that for both 71 = n" and tt = ir', 7c ab (x) — fr ab (x) — > and for both 
g = g" and g = g', g a b( x ) —g a b( x ) — ► 0, all terms vanishing sufficiently fast as 
\x\ — ► 00. In this case g = exp(7 T /2)^exp(7/2). Note that the asymptotic 
fields can depend on x. In this case the coherent-state overlap still holds 
in the form given. This kind of asymptotic behavior reflects a change of 
the fiducial vector \r)), which now depends on the explicitly chosen asymp- 
totic form for the momentum and metric — or, equivalently, as we effectively 
do, one can hold \rj) fixed and change the representations of the operator 
[cf., (71)]. By choosing a suitable asymptotic momentum and metric one 
can, in effect, redefine the topology of the underlying space S. However, for 
simplicity, we shall assume simple Euclidean-like asymptotic behavior of the 
momentum and metric [ n ab (x) = and Y s ( x ) — 0> i- e -> 9ab( x ) = Sab]. 

A second way to study the coherent-state overlap is under coordinate 
transformations. Observe that (iT",g"\7r',g') is invariant if, everywhere, we 
make the replacements 

(i) b(x) o?x by b(x) cfe = b(x) d 3 x, 

(ii) g k \x) by g kl (x) = M k (x) g rs (x) M l s (x), 

(iii) 6(x) _1 7r w (x) by bix)- 1 ^ 1 (x) = b^x)' 1 M k r (x) n rs (x) M l s (x) , 

all for an arbitrary nonsingular matrix M = {M k }, M k (x) = (dx k /dx r )(x), 
that arises from a nonsingular coordinate transformation x — * x = x(x). It 
suffices to restrict attention to those coordinate transformations continuously 
connected to the identity. When S is compact, a wide class of M is allowed; 



when S is noncompact, the allowed elements M must also map coherent 
states into coherent states for the same fiducial vector. This restriction ex- 
cludes any connection by coordinate transformations of two field sets with 
fundamentally different asymptotic behavior; such fields live in disjoint sets. 
The invariance under admissible coordinate transformations is symbolized by 
the statement that 

(W",g"\w',g') = (n",g'y,g f ) (40) 

for all suitable M. Since the allowed M form a representation of the con- 
nected component of the diffeomorphism group, it follows from this identity, 
and suitable continuity, that for sufficiently restricted M, the transformation 
\ir : g) — > |tF, #) is induced by a unitary transformation, specifically that 

\W,g) = U(M)\<K,g), (41) 

U{M) = exp [-ijN j (x) Hj (x) al 3 x] . (42) 

Here T~ij(x) denotes a local operator tensor density of weight one while N^(x) 
denotes a c-number tensor with sufficiently rapid decay at spatial infinity. 
Furthermore, using the shorthand that jN^Hj = jN^(y)Hj(y) cfiy, the con- 
nection between M and iV- 7 is implicitly given by 

M"rt = g ab - {JNJH J} g ab } + ^{jN k H k , {JN^H„ g ab }} + ... 

sfif^-y, (43) 

where {■, ■} denotes the classical Poisson brackets, and specifically, e.g., 
{g ab (x), ir rs (y)} = S r a lS(x,y). In this expression, H,(x) = -£j fe (x)7r*^(x), 
j = 1,2,3, where ( )u is the covariant derivative with respect to the 3x3 
metric, denotes the classical generators of the diffeomorphism group [2"6fl. The 



relationship of Hj(y) and Hj(y) may be determined as follows. Expansion of 
the relation 

{K>\g"\e-^ N3 ^(x)^ W , g >) = (7r",<7V,^> (44) 

to first order in A" leads to 

(7i",g"\fN\x)H,(x)d 3 x\n',g')/(n",g'y,g') 
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-i fb(x) d 3 x([g" kl (x) + g ,k \x)} +ib(x)- 1 [n" k \x) - n' kl (x)]) 
x[8g' kl (x) - ibix)- 1 87r' kl (x)] 
+i \ fb(x) d 3 x g' kl (x) 5g' k \x) , (45) 



where 

8g' kl {x) = g' k ]{x) N j (x) - g' j \x) N k 3 (x) - g' kj {x) N l tj (x) , (46) 

and likewise for 5n' kl (x). This relation determines the coherent state matrix 
elements of 7ij(x). Finally, observe that the diagonal coherent state matrix 
elements read 

(n,g\ Hj{x) \n,g) = -g^^n^x) (47) 

in conformity with the weak correspondence principle. 

3.3 Path integral construction 

If the given coherent states \ir,g) possessed a resolution of unity, namely a 
nonnegative measure p(ir, g) (countably or even finitely additive) such that 

f W , g"\n, g}(n, g\n' , g') dp(n, g) = (it", g"\ir', g') , (48) 

then the construction of a path integral for the reproducing kernel would 
be straightforward and would follow the pattern illustrated in Sec. 2 for 
a single degree of freedom. However, for the proposed reproducing kernel 
(n",g"\iT',g'} given in (39) no such measure exists and thus the traditional 
resolution of unity is unavailable. Consequently, as defined, {\tt, g)} is a set 
of weak coherent states. 

A similar kind of problem arose in the simple model discussed in Sec. 2 
(when < f3 < |). In that case, the construction of a path integral represen- 
tation proceeded in an alternative manner beginning first with a polarization. 
We assert that each of the given Hilbert space representatives, 

N 

^(t, 9) = (n,9\il>) = ^2 a n (tt, g\ir n , g n ) E C , (49) 

n=l 
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satisfies the functional differential equation [cf., (17)] 
Bl{x)^,g) = 



ig rt (x) 



A A 

+ S r s + b(x)" l g st (x) 



Sir ts (x 



Sgtr(x) 



for all spatial points x. Next, let us introduce the operator 
A=\j B s r (xY B r s (x) b(x) d 3 x , 



ip(7r,g) = 
(50) 

(51) 



and observe that A > 0. Thus, with T > and as v — > oo, it follows that 
J{y)e~ vT - A , for some J{y), serves to select out the subspace where (50) is 
fulfilled. Just as in the toy model of Sec. 2, the operator A is a second 
order (functional) differential operator, and, as a consequence, a Feynman- 
Kac-Stratonovich path (i.e., functional) integral representation may be in- 
troduced. In particular, we obtain the formal expression [cf., (21)] 

(n",g"\n',g')= lim AT f exp[-i J g ab ir ab d 3 x dt] 

u^oo J 

x exp{-(l/2i/) |[&(x)-W^ 6c ^ a + K^VWa] d 3 xd£} 

UU^t)^,*). (52) 



X 



x,t a<6 



Here, let us interpret t, < t < T, as coordinate "time". On the right-hand 
side the canonical fields are functions of space and time, that is 



9ab = 9ab{x,t) , 71 



ab 



Ti 



ab , 



X,t) 



(53) 



the overdot ( ' ) denotes a partial derivative with respect to t, and the integra- 
tion is subject to the boundary conditions that n(x, 0), g(x, 0) = tt'(x), g'(x) 
and tt(x, T), g(x, T) = tt"(x), g"{x). Observe that the field tt enters this path 
integral expression only in the form n; this fact is responsible for the result of 
the path integral depending only on n" — it' . It is important to note, for any 
v < oo, that underlying the formal measure given above, there is a genuine, 
countably additive measure on (generalized) functions gu and -n rs . Loosely 
speaking, such functions have Wiener-like behavior with respect to time and 
(5-correlated, generalized Poisson-like behavior with respect to space. 

While (52) is invariant under spatial diffeomorphisms, it is less evident 
that it is also invariant under transformations of the time coordinate (by 
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itself) .0 Formally speaking, the role of the limit v — > oo is to remove the 
effects of the continuous-time regularization. It is clear, however, that there 
is no need that removing those effects must be done in a uniform way in- 
dependent of x. Thus we may replace v by vN(x) — now under the integral 
sign — where N(x), < N(x) < oo, is smooth and reflects the relative rate 
at which the regularization is removed at different spatial points. The end 
result is invariant under such a change. Moreover, at each point x we can 
run the process with different "clock" rates, i.e., N(x) dt — > N(x, t) dt so long 
as the elapsed time is qualitatively unchanged. This remark means that we 
can choose any smooth lapse function N(x,t), < N(x,t) < oo, with the 
consequence that 



<y 



,g"\n',g') = lim M' / exp[-i f g ab 7r ab d 3 x dt} 

u^co J ■ 

x ^ V {-(l/2v)$[b(x)- l g ah g cd K bc it da + b(x)g ab g cd g bc g da ] N{x,t)- X d 3 xdt} 

lH[dn ab (x,t)dg ab (x,t). (54) 



x 

x,t a<b 



The necessary conditions for this more general expression to hold are, for all 
T, < T < oo, and at all x, that 

fl N(x, t) dt < oo , (55) 

f™N(x,t)dt = oo. (56) 

In this sense we observe that our formal path integral representation (52) for 
the coherent-state overlap is actually invariant under transformations of the 
time coordinate. 



3.4 Metrical quantization 

The formal, ^-dependent, weighting factor in the path integral expression 
(52) involves a metric dT? on the classical phase space which may be read 
out of the expression 

dY?/dt 2 = Jib' 1 g kl g rs 7r lr 7r sk + bg kl g rs g lr g sk ] d\ . (57) 



The author thanks A. Ashtekar for raising the question of temporal transformation 
properties. 
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As presented, this expression for dYJ 2 is a derived quantity. Alternatively, it 
is clear that one could start the analysis by postulating a specific functional 
form for cffi 2 to be used in a continuous-time regularization in the path inte- 
gral construction of the reproducing kernel, and finally, by appealing to the 
GNS (Gel'fand, Naimark, Segal) Theorem ||27|| , to recover the representation 
of the local field operators o~ki{x) and K r s (x). Adopting a metric on the classi- 
cal phase space as the first step in a quantization procedure is called metrical 
quantization |28|] . To carry out such a scheme for gravity, it is necessary that 
any postulated gE 2 satisfy several properties. First, it must be diffeomor- 
phism invariant and second, on physical grounds, it should only depend on 
dir kl (or -k kl for dY? /dt 2 ) and not on tt m itself. Hence, we are initially led to 
consider 

dY 2 /dt 2 = J[b- 2 L abcd 7r bc 7r da + M abcd g bc g da ] b(x) d 3 x (58) 

for suitable, positive-definite tensors L and M constructed just from g ki . 
The given choice for L and M, i.e., L abcd = \\g ab g c d + g ac gbd] and M abcd = 
\[g ab g cd + g ac g bd ] satisfy M = L~ x as matrices. This choice is very natural and 
moreover is identical to the form suggested by the study of certain GL(3,R) 
coherent states for a 3 x 3 positive-definite matrix degree of freedom 



Nevertheless, in a metric-first quantization scheme, it is appropriate to 
examine other choices as well. For example, a term such as b(x)~ l g ab (x) , k ab (x) 
might be included, but this term may be eliminated by a translation of the 
momentum. Additionally, one may consider nonlocal contributions involving, 
for example, the term g bc (x)gd a (y), together with a kernel K(x,y) specifying 
the interrelationship of the field at y to the field at x. However, no satisfactory 
solution for K(x, y) other than one proportional to 5(x, y) will lead to an 
expression for dY? that is invariant under all diffeomorphisms. In point of 
fact, the possible choices for L and M are rather limited, especially when 
one requires that the (formal) integration measure at each point is canonical 
and thus has the form n a <;, d7r ab dg ab . For example, for A > 0, let us consider 
the proposal that 

Labcd (A) = \[g a b9cd + 9ac9bd + (A - l)g bc g da ] . (59) 

Then in order to lead to a canonical integration measure it would be necessary 
that 

M abcd (X) = \{g ab g cd + g ac g bd + (A" 1 - l)g bc g da ] . (60) 
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Only for A = 1 is M = L which is just the choice we have made. (The form 



for the DeWitt metric |29|], where A = — 1, is excluded because we require 
that L and M be positive definite.) 

The preceeding discussion has rather convincingly suggested the specif- 
ically chosen functional form for dT? — apart from one issue. It may seem 



even more natural to choose 30 



dY?/dt 2 = j[g- 1 / 2 g ab g cd 'k bc 7c da + g 1/2 g ab g cd g bc g da } d 3 x (61) 

rather than the choice we have made. This is a natural choice from a classical 
point of view, but it is less satisfactory from a quantum point of view. In 
either case, observe that a path integral such as (52) involves fields with 3 + 1 
independent variables; however, there are no space derivatives involved, only 
time derivatives. Such a model is known as an ultralocal quantum field theory, 
and by now there is much that is known about the rigorous construction and 
evaluation of such nontrivial (i.e., non-Gaussian) functional integrals through 
the study, for example, of ultralocal scalar quantum fields |3T|. It is through 



the analysis of the gravitational models as ultralocal quantum field theories 
that the metric (61) is ruled out; for a simple reason, see Section 5. 

3.5 Operator realization 

In order to realize the metric and scale fields as quantum operators in a 
Hilbert space, Sj, it is expedient to introduce a set of conventional local 
annihilation and creation operators, A(x, k) and A(x, k)\ respectively, with 
the only nonvanishing commutator given by 

[A{x, k), A{x', k'Y] = 5{x, x') S(k, k')l , (62) 

where 1 denotes the unit operator. Here, x G R 3 , while k = {k rs } denotes a 
positive-definite, 3x3 matrix degree of freedom confined to the domain where 
{k rs } > 0. We introduce a "no-particle" state |0) such that A(x, k) |0) = for 
all arguments. Additional states are determined by suitably smeared linear 
combinations of 

A(x u hY A(x 2 , k 2 y ■ ■ ■ A(x p , k p y |0> (63) 

for all p > 1, and the span of all such states is Sj provided, apart from constant 
multiples, that |0) is the only state annihilated by all the A operators. Thus 
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we are led to a conventional Fock representation for the A and A' operators. 
Note that the Fock operators are irreducible, and thus all operators acting 
in f) are given as suitable functions of them. 

Next, let c(x, k) be a possibly complex, c-number function and introduce 
the translated Fock operators 

B(x,k) = A(x,k) + c(x,k)l, (64) 

B(x,kY = A(x.ky + c(x,k)*l. ( 65 ) 

Evidently, the only nonvanishing commutator of the B and B^ operators is 

[B(x, k), B(x', k'Y] = S(x, x') 5{k, k')l , (66) 

the same as the A and A* operators. With regard to transformations of the 
coordinate x, it is clear that c(x, k) (just like the local operators A and B) 
should transform as a scalar density of weight one-half. Thus we set 

c{x,k) = b{x) 1 / 2 d{x,k) , (67) 

where d(x, k) transforms as a scalar. The criteria for acceptable d(x, k) are, 
for each x, that 

J + \d(x,k)\ 2 dk = oo , (68) 

f + k rs \d(x,k)\ 2 dk = 25 rs , (69) 

the latter assuming [cf., the discussion following (39)] that (jki(x) = Ski- 
In (68) and (69) we have introduced dk = IT a < fe dk a b, and the symbol "+" 
signifies an integration over only those k values for which {k^} > 0. 
We shall focus on only one particular choice for d, specifically, 

J£ e -tr(fc) 

dix ' k ^iMW (70) 

which is everywhere independent of x; K denotes a positive constant to 
be fixed later. The given choice for d corresponds to the case where the 
asymptotic fields tx m {x) = and (jki(x) = Ski- [Remark: For different 
choices of asymptotic fields it suffices to choose 

d(x,k) ->d(x,k) = -j— — . (71) 
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We shall not explicitly discuss this case further.] 

In terms of these quantities, the local metric operator is defined by 

a ab (x) = b{x)- 1 j + B(x, ky k ab B{x, k) dk , (72) 

and the local scale operator is defined by 



kKx) = -i±f B(x, kY {k st ~£) tr - IF*!*,) B(x, k) dk . (73) 



s » 



Here d = d/dk st , d = d/dk rt acting to the left, and a ab (x) transforms 
as a tensor while K r s (x) transforms as a tensor density of weight one. It 
is straightforward to show that these operators satisfy the required affine 



commutation relations, and moreover that [31, 32 



(01 e i J ^^abW'Px e -ij ^(x)K, r s (x)d 3 X |Q\ 

= exp{-K 2 fb{x) d\j[e~ 2&abKb - e -^^)WK*) e -[(^+s a6 (*))M] d k/{det k) 2 } 
= exp[-2/6(x) d 3 xln{[det{g ab {x))} 1/2 det{±[5 ab + g a \x)] - i\b(x)- l n ab (x)})} , 

(74) 

where K has been chosen so that 

K 2 f + k rs e" 2tr(fc) dife/(det k) 2 = 25 rs . (75) 

An obvious extension of this calculation leads to (39). 

3.6 Local operator products 

Basically, local products for the gravitational field operators follow the pat- 
tern for other ultralocal quantum field theories |3l], [32] . As motivation, con- 
sider the product 

o- a b(x)a cd (y) 

= b(x)~ 2 JJ + B(x, k)^ k ab [ B(x, k), B(y, k'y]k' cd B(y, k') dk dk' 

+ : a ab {x)a cd (y) : 
= b(x)~ 2 5(x, y)J + B(x, k) ] k ab k cd B(x, k) dk+ : a ab (x)a cd {y) : , (76) 
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where : : denotes normal ordering with respect to A and A*. When y = x, 
this relation formally becomes 

cr a bi.x)a cd (x) = b(x)~ 2 5(x,x)f + B(x,kyk ab k cd B(x,k)dk+ : a ab (x)a cd (x) :(77) 

We define the renormalized (subscript U R") local product 

[cr a b{x)<J cd (x)} R = b(xy 1 j + B(x, kY k ab k cd B(x, k) dk (78) 

after formally dividing both sides by the divergent dimensionless "scalar" 
b(x)~ l 5(x, x).0 Higher-order local products exist as well, for example, 



<T aibl (x)cr a2 2 (x)a a3b3 (x) ■ ■■a apbp (x)] 



h 



= b(x)- l j + B(x, k)^ {k aibl k a ^k a:ib:i ■ ■ ■ k apbp ) B(x, k) dk , (79) 

which, after contracting on b\ and b 2 , implies that 

[cr aib (x)a a2 \x)a azb3 (x) ■ ■ ■ cr apbp (x)] R = 5* 2 [cr asba (x) ■ ■ ■ a apbp (x)] R . (80) 

It is in this sense that [cr a b(x)a bc (x)} R = 5 a . 

We take up only one further point regarding local products. It is rather 
natural |y], HU to try to define the local momentum "operator" by 



» rs (y) = -i\J + B{y, ky (T b - -T b ) b( v , k) dk , (si) 

but this expression only leads to a form and not a local operator. Further- 
more, the putative canonical commutation relation becomes 

Mx), » rs (y)] = 1 52 6(x, y) b(x)- 1 J + B(x, k) B(x, k) dk 

= 1 52 5(x,y)[f + \d(x,k)\ 2 dk + ...], (82) 

which has a divergent multiplier and is, therefore, not even a form. On the 
other hand, it is true that 



l[a r i(x)fi ls (x) + fi sl (x)a lr (x)} R 

yii 

= -i±f + B(x,k) (k rl d ' 

= 4(x) (83) 



lis i s i 

\J ,B(x, k) (k r i d d ki r ) B(x,k)dk 



as claimed. 



5 For scalar ultralocal theories, the formal dividing factor is the divergent dimensionless 
"number" b~ 1 6(0), where b > is an arbitrary factor with suitable dimensions. For gravity, 
b — > b(x), our scalar density of weight one. Note that limits involving test functions offer 
a rigorous definition of the renormalized product plL B3] . 
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4 Imposition of Constraints 

Gravity has four constraints at every point x G S, and, when expressed 
in suitable units, they are the familiar spatial and temporal constraints, all 
densities of weight one, given by |26 



'a{x) -- 


= -g ab (x)n bc lc (x) , (84) 


'(*) = 


= \g{xy 1/2 [g ab (x)g cd (x) + g a d(x)g cb (x) - 2g ac (x)g bd (x)} 




X7i ac (x)n bd (x) + g(x) 1/2 {3) R(x) . (85) 



The spatial constraints are comparatively easy to incorporate since their gen- 
erators serve as generators of the diffeomorphism group acting on functions 
of the canonical variables. Stated otherwise, finite spatial diffeomorphism 
transformations map any coherent state onto another coherent state as in 
(41) and (42). However, this is decidely not the case for the temporal con- 
straint. What follows is an account of what to do about these constraints in 
principle; in Part II on this subject, we will discuss how to accomplish these 
goals. 

One satisfactory procedure to incorporate all the necessary constraints 
is as follows. Let {h p (x)}^ =1 denote a complete, orthonormal set of real 
functions on S relative to the weight b{x). In particular, we suppose that 

Jh p (x) h n (x) b(x) d\ = 5 pn , (86) 

b(x)E™ 1 h p (x)h p (y) = 6(x,y). (87) 

Based on this orthonormal set of functions, we next introduce four infinite 
sequences of constraints 

H (v)a = Jh P (x) H a {x) d 3 x , (88) 

H {p) = Jh p (x) H(x) d 3 x , (89) 

1 < p < oo, all of which vanish in the classical theory. 

For the quantum theory let us assume, for each p, that 7i( p ) a and TCm 
are self adjoint, and even stronger that 

Xl = S^ =1 2-p [Zl =1 (H {p)a f + (H {p) ) 2 ] (90) 
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is self adjoint for all P < oo. Note well, as one potential example, the 
factor 2 _p introduced as part of a regulator as P — > oo; we comment on this 
regulator in the next section. For each 5 = 5(h) > 0, let 

E P = E(Xp < 5 2 ) (91) 

denote a projection operator depending on Xp and S as indicated. How such 
projection operators may be constructed is discussed in |7| and will be dealt 
with in Part II. Let 

S P = \im sup (7!-,g\E P \iT,g) , (92) 

which satisfies Sp > since Ep ^ when restricted to sufficiently large S. 
Finally, we define 

((<*", g>y,g')) = limsup Sp 1 (n",g"\E P \n',g') (93) 

as a reduction of the original reproducing kernel. The result is either trivial, 
say if S is too small, or it leads to a continuous, positive-definite functional 
on the original phase space variables. We focus on the latter case. 

To obtain the final physical Hilbert space, one must study ((it", g"\ir', g')) 
as a function of the regularization parameter 5. Since gravity has an anomaly 



33|], there should be a minimum value of 6, which is still positive, that defines 
the proper theory, rather like the example in (30). Assuming we can find and 
then use that value, ({it", g"\^' ■, g')) becomes the reproducing kernel for the 
physical Hilbert space i3 p h y s- Attaining this goal would then permit the real 
work of extracting the physics to begin. 

Our discussion regarding constraints in this paper has indeed been brief. 
Although the program we have in mind is not simple, it has the virtue of 
being realizable, at least in principle. After all, before any calculational 
scheme is developed it is always wise to ensure that the object under study 
has a good chance of existing! 

5 Discussion 

In the preceding sections we have outlined an approach to quantum gravity 
that it is somewhat different than currently considered. As background to 
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our philosophy, let us briefly review some of the common weak points in the 
standard ways of quantizing gravity and use these comments as motivation 
for our approach. 

5.1 Traditional viewpoints & commentary 

(i) Viewed by way of conventional perturbation theory, quantum gravity has 
two main difficulties of principle. On the one hand, the perturbative split 
of the metric in the form g^ u (x) = r)^ + h^ v (x) (or any other background 
metric), with canonical quantization of the "small" deviation h^ u (x) violates 
signature properties since in that case the spectrum of h^ v (x) is unbounded 
above and below. On the other hand, as an asymptotically nonfree theory, 
gravity is nonrenormalizable and poorly described by a perturbation theory 
which needs an unending addition of distinct counterterms with divergent 
coefficients. 

To address these obstacles, we first note that the affine approach guar- 
entees a proper metric signature from the very beginning, and second we 
remind the reader that certain asymptotically nonfree, nonrenormalizable 



models have indeed been solved |3T| , and their solution procedures form the 
core of the present approach to quantize gravity. 

(ii) While the constraints of classical gravity are first class, there is an 
anomaly in the quantum constraints and thus they are effectively second 
class. Usual views toward second-class constraints involve solving and elim- 
inating them, introducing and then quantizing Dirac brackets, or the con- 
version of second-class constraints into first-class constraints. Each of these 
methods is often complicated and not all are guaranteed to be valid be- 
yond a semiclassical treatment if the classical constraint hypersurface has 



a non-Euclidean geometry [pl |. These difficulties have stimulated searches 



to get around the second-class character altogether, either by introducing 
non-Hermitian constraint operators that may close algebraically p5[ , or by 
introducing additional fields and space-time dimensions until the anomaly 
cancels. 

Regarding these comments, we accept the anomaly and the second-class 
constraints that it implies. Giving up a classical symmetry is is not so hereti- 
cal as it may seem. For example, Hamiltonian classical mechanics enjoys a 
full covariance under general canonical coordinate transformations, but that 
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invariance is not preserved in its classical form when we go to the quan- 
tum theory. For example, consider the classical Poisson brackets for a set of 
generator elements 

{e ap+bq , e cp+dq } = (be - ad) e ( a + c )p+(b+d) q ^ ^ 

where a, b, c, and d are parameters, while in canonical quantum mechanics 
we have the corresponding commutator algebra 

[e aP+bQ , e cP+dQ ] = (2i) sin[(6c - ad)/2] e i^)p+(b+d)Q _ (Q5) 

These expressions agree in their algebraic structure for selected elements, but 
not for the whole algebra. Equivalence begins to break down at the quadratic 
level, which is exactly the case for the temporal constraint in gravity. In par- 
ticle mechanics there are sound physical reasons []36] for this "breakdown" of 
symmetry, and attempts to restore the symmetry — as in geometric quantiza- 
tion pTj — go counter to such sound physical principles. There is no reason 



that a similar scenario does not hold for gravity. The breakdown of the clas- 
sical symmetry and the appearance of a quantum "anomaly" (better called a 
"quantum mechanical symmetry breaking" [J38|1 ) could, just as in the quan- 
tum mechanics case, carry real physics. 

Accepting the second-class nature of (part of) the constraints of quantum 
gravity means a different approach must be taken. As already noted, earlier 
approaches required solving for and eliminating the unphysical variables, the 
introduction of Dirac brackets, etc., all of which are rather technical and may 
be extremely complicated. In the present view, afforded by the projection 
operator approach to constraints — second-class constraints included — none of 
these particular complications arise. Instead, one projects onto the state (or, 
with degeneracy, states) for which the sum of the square of the constraints is 
bounded. Why the square and not the fourth power? Using the fourth power 
would not be wrong; the only change would involve a unitary transformation 
of the original result, which maps one set of "ground" states onto another set 
of "ground" states. The square is chosen for simplicity, not for any reasons 
of exclusivity. 

5.2 Relation to earlier work 

Pilati, in a series of papers || (see also ||), analyzed a strong coupling 
model of quantum gravity in which the temporal constraint given in (85) 
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was modified to read 

H'(x) = \g(x)- 1,2 [g ab (x)g cd (x) + g a d(x)g cb (x) - 2g ac (x)g bd (x)] -K ac (x)ii bd (x) , 

(96) 

namely, the second term involving the scalar curvature ^R(x) based on the 
metric g a b{x) was dropped. The reason for doing so was to achieve a theory in 
which the temporal constraint H'(x) itself was patterned after the Hamilto- 
nian density of an ultralocal theory. This modification was thought to be ad- 
vantageous because then all the machinery developed for ultralocal quantum 
field theory could be used for the strong coupling gravitational model. Once 
that model was under control, it was the hope to reintroduce the dropped 
term by a perturbation theory analysis. Unfortunately, the reintroduction of 
dropped terms involving spatial gradients has never been successfully accom- 
plished by a perturbation analysis about a non-Gaussian ultralocal model. 
This failure is most likely because such "interaction terms" generally amount 
to nonrenormalizable perturbations of the unperturbed (ultralocal) models. 

The program advanced in the present paper takes a different view toward 
these issues. 

First, we focus on kinematics with the knowledge that for pure grav- 
ity the Hamiltonian operator vanishes, as it does in any situation that is 
reparametrization invariant. In its place we find constraints, and the real 
physical content of the theory lies in the particular constraints. However, 
before the constraints can be introduced, there must be a "primary con- 
tainer" to receive them. In our case, this primary container is the Hilbert 
space and set of relevant operators prior to the introduction of any form of 
the constraints, and which is based on the fundamental physical nature of 
the variables, i.e., positive-definite, 3x3 matrix- valued, local field operators, 
etc. This is the preferred procedure: Quantization before the introduction 
of any constraints. At this primary level, there is no coupling of one degree 
of freedom with another — any coupling comes through the enforcement of 
specific constraints. Hence, in the primary container the degrees of freedom 
are mutually independent of each other. For finitely many kinematical de- 
grees of freedom this means that the Hilbert space is a product over spaces 
for each of the separate degrees of freedom; in a field theory, this indepen- 
dence means that the kinematical operators enter as ultralocal field operators. 
Consequently, even though the several constraint operators waiting to be in- 
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troduced may themselves not be ultralocal in nature, the primary container 
itself, which has been prepared to receive them, is ultralocal. 

At this point the reader may wish to reexamine (39) — in essence our 
"primary container" — to recall the appearance of an ultralocal state on a set 
of field operators. Apart from the 3x3 matrix character, the functional 
form of (39) emerges from (i) the product of iV expressions of the form (16) 
for independent arguments p n ,q n ,r n ,s n , 1 < n < N, (ii) the replacement 
of j3 by 6 n A and (p n — r n ) by (p n — r n )A, and (iii) the limit as N — > oo, 
A — > such that £(•) b n A^ f (•) b(x) dx. In this way we have created, from 
a collection of independent single affine degrees of freedom, the reproducing 
kernel for affine gravity in 1 + 1 dimensional space. In a similar manner, a 
set of independent 3x3 affine degrees of freedom can be (and were) used to 
build an ultralocal representation for 3x3 metric and momentum fields in 
(39); moreover, this type of construction does not favor the "natural" phase- 
space metric (61). In summary, we emphasize that whenever the "dynamics" 
appears through constraints, the primary container should be ultralocal in 
character. We next turn our attention to the introduction of the constraints. 

In the projection operator approach it is recognized from the outset that 
the physical Hilbert space — or better the regularized physical Hilbert space — 
is a subspace of the original Hilbert space that is uniquely determined by an 
associated projection operator E. Whatever form the constraints may take, 
they are "encoded" into the projection operator E, and a regularization 
means that the constraints are satisfied to a certain level of precision deter- 
mined by a regularization parameter 5. How to turn constraint operators 
into projection operators in general has been discussed in 0; as regards the 
gravitational case, that project will be discussed in Part II . 

Continuing still in a general framework, let us consider an expression that 
may be used either to generate dynamics or to enforce constraints. From a 
classical point of view, and especially from a path integral point of view, it 
may seem that quantities used in either of these ways may be rather similar. 
However, it is important to already understand that there is a fundamental 
distinction between the use of a quantum operator either (i) to generate 
unitary transformations or (ii) to serve as a constraint operator in a given 
system. In the first case, the operator must be self adjoint and thus densely 
defined, while in the second case the operator may be defined on only the zero 
vector! This fact has profound consequences. In particular, to have a self- 
adjoint generator requires that the operator representation in the primary 
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container must already be finely "tuned" to ensure that the generator which 
will be introduced is self adjoint (as in Haag's Theorem |]39|). For constraint 
operator imposition this need not be the case, and the reason this is so 
is because we allow for changes, i.e., adaptions, of the primary container 
representations through the process of reduction of the reproducing kernel. 
As an example, let us consider only the local temporal operator 7i(x) for 
gravity. On the one hand, to generate unitary time evolutions it may be 
necessary that f7~C{x) d?x be self adjoint. On the other hand, to enforce 
constraints, it is only necessary that H( p ) [cf., (89)] be "small", but there 
is no requirement that these operators must be uniformly "small". Instead 
they can be "small" in the sense that X 2 = E^ =1 s p (7i( p )) 2 is "small", where 
the set of positive constants {s p } serve as regulators to control convergence 
of the series. The example chosen was that s p = 2~ p , but there is nothing 
special about that choice. Any reasonable choice that leads to a self adjoint 
operator X should lead to the same reproducing kernel in the final analysis 
when the regularization parameter 5 attains its final value for the problem 
at hand. 

A rather simple example of the general procedure discussed above can 
be seen in studies involving product representations jE|. Additionally, it is 
instructive to reanalyze the relativistic free field by this procedure to see how 
ultralocal representations turn into nonultralocal representations. 

Suffice it to say, it is this vast difference between the required nature of 
constraint operators and unitary generators that permits us to start with 
ultralocal field operator representations and emerge with nonultralocal op- 
erators in the physical Hilbert space. 
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